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lent atmosphere. 


Unci  ass  1 f1c<l 


IMAGING  OF  JNCOHERENT^GBJECTS  BY  J 

[OTBNSITS  CORRELATION  MTERPEROMETRY^,! 


Wwo“a'?oTLSr‘  ■’<’»'■  •"« 

Ei=.ribu.i::  s ^ <’ 

A.  D,  BLOSB 

rochalcal  tefonMtion  omcer 


/f  j 


period:  1 Oot.-'l971  - 30  Sept^_1976  ^ 

principal  Investigato^j^r.^H^^am_/ 

.:.  J/y  Electrical  Engineering 


■P:2  ! 


TABLE  OF  CONTENTS 


1.  Introduction 

2.  Laboratory  Simulation 

3.  Field  Experiment 

4.  Computer  Simulation 

5.  Conclusion 
Reference 

I Appendix:  A.  A general  purpose  interactive  experimental 

t *■ 

I control  system 

I 

I « 

I B.  Review  of  pertinent  research 


i V 


1.  INTRODUCTION 


The  objective  ox  this  research  was  to  develop  a new  electro- 
optical  imaging  technique  for  incoherent  light  source,  which  is 
not  seriously  disturbed  by  the  atmospheric  turbulence. 

The  basic  physical  quantity  to  be  processed  for  imaging  of 

incoherent  light  sources  is  the  mutual  intensity  as  a function  of 

distance  between  two  observation  points.  The  mutual  intensity 

is  the  crcss-correlation  of  incident  complex  wave  amplitudes 

transmitted  from  an  Incoherent  light  source.  The  tine  delay 

between  these  two  wave  amplitudes  are  taken  to  be  zero,  where  the 

magnitude  of  cross-correlation  is  maximum.  The  mutual  intensity 

of  an  Incoherent  light  source  is  generally  a complex  quantity, 

specified  by  both  magnitude  and  phase.  When  the  mutual  intensity 

is  determined  as  a function  of  the  distance  between  two  observation 

points,  the  angular  intensity  distribution  of  the  incoherent  object 

can  be  obtained  from  the  mutual  intensity,  based  on  the  Van  Cittert- 
2 

Zemike  theorem.  In  the  far  field  case,  the  angular  Intensity 
distribution  is  given  by  the  Fourier  transform  of  the  mutual  intensity. 

The  usual  photographic  image  of  an  incoherent  source  by  using 
a lens  system  can  be  interpreted  as  the  results  of  analog  Fourier 

3 

transform  of  mutual  intensity  on  the  pupil  plane. 

. Under  the  presence  of  the  atmospheric  turbulence,  the  photo- 
graphic image  of  an  incoherent  object  will  be  seriously  degraded 
by  the  phase  and  amplitude  fluctuations  as  well  as  the  fluctuation 
of  angle  cf  Incidence.  Thus,  the  method  of  indirect  imaging  of 
an  incoherent  object  through  turbulent  atmosphere  based  on  the 
intensity  correlation  interferometry  and  image  processing  becomes 
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an  attractive  alternative,  worth  further  investigation.  Both  the 

intensity  correlation  interferometer  with  a coherent  background 

and  the  triple  intensity  correlation  interferometer  proposed  by 

the  principal  investigator  can  determine  both  magnitude  and  phase 

of  the  mutual  intensity,  without  being  seriously  disturbed  by  the 

4 5 

atmospheric  turbulence.  ' It  should  be  noted  that  the  original 
Hanbury  Brown-Twlss  intensity  interferometer  is  not  sufficient  for 
general  purpose  imaging  because  the  phase  of  mutual  intensity 
cannot  be  determined  by  the  interferometer. 

The  essential  feature  of  the  indirect  imaging  based  on  the 
modified  intensity  correlation  interferometer  is  that  the  mutual 
intensity  can  be  determined  without  being  seriously  disturbed  by 
the  atmospheric  turbulence. 

This  type  of  imaging  system  is  insensitive  to  the  phase 
fluctuations.  The  amplitude  fluctuations  due  to  atmospheric 
turbulence  will  make  the  signal-to-nolse  ratio  of  detected  mutual 
intensity  smaller,  and  consequently  the  imaging  process  requires 
longer  observation  time  than  the  usual  photographic  imaging. 

After  applying  the  signal  averaging  procedure  to  the  mutual 
intensity  measured,  we  should  be  able  to  recover  the  image  of  an 
original  incoherent  source,  irrespective  of  the  presence  of 
atmospheric  turbulence. 

The  development  of  imaging  system  based  on  the  generalized 
intensity  correlation  interferometer  is  a challenging  research 
topic  in  the  fundamental  optics  as  well  as  in  optical  electronics, 
because  it  requires  new  types  of  opto-electronic  detector  system 
and  real  time  digital  signal  processing. 
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First,  we  report  the  laboratory  simulation  experiments  using 
a pseudo-thermal  source,  the  laboratory  generated,  homogeneously 
turbulent  atmosphere  and  a linear  array  of  photodetectuis.  Second, 
the  field  experiments  of  optical  transmission  through  the  turbulent 
atmosphere  using  a reflector  on  the  campus.  Third,  some  results  of 
preparing  computer  simulation  of  optical  transmission  through  tur- 
bulence are  described,  with  emphasis  on  the  generation  of  non- 
gaussian  random  processes  with  prescribed  moments  or  cumulants. 
Brief  review  of  pertinent  research  works,  and  the  interactive 
experimental  control  system  program  based  on  the  FOCAL  language 
for  FDP  15/20  are  described  in  the  appendix. 
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2.  LABOIIATORY  SIMULATION 


/k  preliminary  experiment  on  intensity  correlation  interfero- 
metry was  performed  using  a pseudo- thermal  source  as  an  incoherent 
object  and  analog  multiplication.  The  pseudo-thermal  source 

O 

consists  of  a rotating  ground  glass  with  the  He-Ne  6328A  gas  laser. 
The  statistics  of  laser  radiation  scattered  by  a rotating  ground 
glass  were  measured  by  using  the  experimental  system  illustrated 
in  Fig.  1,  in  which  an  Analog-to-Digital  convertor  interfaced  to 
the  PDF  15/20  digital  computer  played  the  important  role.  The 
histogram  of  irradiance  fluctuations  obtained  is  very  close  to 
thr'.  exponential  distribution  (cf.  Fig.  2),  This  Indicated  that 
the  isticc  of  wave  amplitude  fluctuations  are  close  to 
I'snes^ian  and  the  pseudo-thermal  source  can  be  treated  as  an  inco- 
herent light  source.  The  schematic  diagram  of  an  intensity 
correlation  interferometer  tested  by  using  an  analog  multiplier 
(Intronlcs  501)  Is  illustrated  in  Fig,  3,  and  a typical  result 
is  shown  in  Fig,  A,^  Since  we  intended  to  develop  the  digital 
signal  processing  system  for  studying  various  intensity  correlation 
interferometers,  the  system  using  intensity  correlation  by  analog 
multiplication  was  not  elaborated. 

Prior  to  developing  the  pseudo-thermal  source  using  a 
rotating  ground  glass,  we  also  investigated  the  statistical 
properties  of  the  light  scattered  by  suspended  Brownian  particles. 
Polystyrene  latex  of  known  particle  diameter  such  as  0,714  ym 
(Dow  Chemical)  were  suspended  in  distilled  water.  The  advantage 
of  this  type  of  source  is  that  it  does  not  require  the  external 
driving  force  such  as  a motor  for  rotating  a ground  glass  plate. 
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Disadvantage  is  that  the  coagulation  of  colloid  may  change  the  long 
range  characteristics  of  the  pseudo-thermal  source.  Hie  dynamic 
scattering  in  nematic  liquid- crystal  and  an  array  of  light  emitting 
diodes  would  also  be  useful  as  a pseudo-thermal  source  for  repre- 
senting the  incoherent  objects. 

Statistics  of  the  pseudo-thermal  source  were  thoroughly 
investigated  by  measuring  moments,  central  moments  and  cumulants 
using  our  digital  instrumentation.  The  details  are  described  in 
the  Ph.D  dissertation  by  Donald  G.  Lubnau  (June  1976).^  We  also 
presented  a paper  on  depolarization  properties  of  random  rough 
dielectric  surfaces  at  the  Optical  Society  of  America  Meeting 
(October  1972),^ 

In  conjunction  with  this  research,  the  principal  investigator 
has  re-examined  the  derivation  of  the  Van  Cittert  theory  of 
partially  coherent  radiation.  By  establishing  the  Fokker-Plank 
equation  for  joint  probability  density  of  partially  coherent 
wave  field,  the  missing  link  between  the  Zemike  theory  and  the 
Van  Cittert  theory  were  clarified.  This  was  presented  at  the 

g 

Annual  Meeting  of  the  Optical  Society  of  America  (October  1972) . 

While  performing  the  measurements  of  irradiance  fluctuations 
by  a photovoltaic  detector  (Si  PN  junction)  we  noticed  that  the 
noise  level  depends  on  the  reverse  bias  voltage  and  intensity  of 
incident  laser  radiation.  Thus,  we  measured  the  excess  noise  in 
silicon  photodiode  (EG  & G,  SGD  100)  while  changing  the  inten- 
sities of  the  Incident  He-Ne  visible  gas  laser  radiation.  This 
is  described  in  the  M.S,  dissertation  by  Rose  A,  Shuttleworth, 
(1974).^ 
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After  several  measurements  of  optical  transmission  through  real 
turbulent  atmosphere,  we  found  It  rather  difficult  to  obtain 
consistent  results  due  to  the  uncertainty  of  atmospheric  conditions. 
Therefore,  we  decided  to  build  a simple  laboratory  generation  of 
the  atmospheric  turbulence  which  will  provide  the  reproductible 
homogeneous,  isotropic  turbulence  over  a reasonably  large  region. 
Consequently,  we  can  now  conveniently  perform  optical  transmission 
measurements  using  the  turbulence  chamber.  This  is  one  of  the 
highlights  of  our  research  activities  during  the  five  years  and 
its  details  will  be  described  in  Arun  K.  Majumdar’s  Ph.D  disser- 
tation (May  1977).^*^ 

The  turbulence  chamber  (100  inch  x 31  Inch  x 9 inch)  consists 
of  10  small  electric  heater/blowers  with  provision  for  heating 
at  750  watts  or  operating  just  fan  only  without  heating  (see 
Figs.  5 and  6).  An  aluminum  foil  screen  and  three  screens  of 
2 rm  altnoinum  wire  meshes  were  placed  in  front  of  the  heaters,  in 
order  to  generate  the  homogeneous,  isotropic  turbulence.  The 
turbulence  within  the  region  (12  inch  x 100  inch  x 3 inch)  at  the 
4 3/4  inch  standard  height  for  optical  measurement  was  found 
locally  homogeneous  and  Isotropic.  The  average  wind  velocity 
and  temperature  were  0,41  m/sec  and  53®  C,  respectively. 

The  turbulence  chamber  was  characterized  by  both  thermal 

and  optical  methods.  The  temperature  structure  function 

2 

D^(r)  * <lT(rj^  + r)-T(rj^)|  > were  measured  in  various  locations 
by  using  the  differential  micro-thermocouple  system  (Chromel/ 
Constantan,  Type  E,  Hy-Cal  Engineering),  (See  Fig.  7,  8 A and  8 B), 
From  the  measurements  of  the  temperature  structure  function  versus 


i 
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probe  distance  r,  we  found  that  the  following  Kolraorgov-Obukltov's 
2/3  power  law  holds  at  probe  distance  from  2 to  6 cm; 


D^(r)  = r^^^ 


The  parameter  for  the  strength  of  temperature  fluctuations 
was  52.9  (®K)^  meter-2/3. 

The  power  spectra  of  temperature  fluctuations  were  obtained 

by  applying  the  Fast  Fourier  Transform  method  to  the  digital 

temperature  fluctuations  stored  in  the  data  acquisition  system. 

(See  Fig.  10).  Within  the  frequency  region  higher  than  5 Hz  but 
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lower  than  80  Hz,  the  theoretical  frequency  dependanace  f of 
power  spectra  for  the  homogeneous  isotropic  turbulence  in  sub- 
inertial  range  is  in  close  agreement  with  our  experimental 
results.  For  an  average  wind  velocity  C.Al  m/sec  we  obtain  the 
inner  and  outer  scales  of  turbulence,  * 5 mm  and  L®  * 8.2  cm, 
respectively. 

The  refractive  index  structure  constants  C were  calculated 

n 

from  measured  by  using  the  formula 


C 

n 


77.6  X 10 


.-6 


0.00753 


P [+ 


,2  ■ ^ S 


2 ' ■ 

The  distribution  of  is  illustrated  in  Fig.  11,  where  the  x 
and  y coordinates  are  parallel  to  the  optical  path  and  airflow, 
respectively. 
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The  strength  of  this  laboratory  generated  turbulence  is 
3 

approximately  10  times  larger  than  that  of  the  real  atmosphere, 
Hiis  is  an  important  feature  of  our  .turbulence  chamber  to  make 
the  simulation  of  the  real  atmosphere  within  the  laboratory 
possible. 

Irradlance  fluctuations  of  a collimated  single  mode  He-Ne 

O 

6328A  laser  beam  were  measured  by  transmitting  the  1"  diameter 
collimated  beam  over  the  2,5  meter  distance  of  the  turbulent 
atmosphere  mentioned  above.  From  the  measured  moments  up  to 
the  fourth  order  of  irradiance  fluctuations  we  obtained  the 
variance  of  log-irradiance  fluctuations  by  using  the  formula: 


C,  (0) 


- <[  In  (I/Io)  - < In  (I/Io)  >1  > 
4 


1 A - 11  , 2 

ss  — ( 6 <y  > -4  <y"^>  + **6  <y> 


1 2 

2 2 2 3 

— — <y  > + 4 <y>  <y  > ••  — <y>  <y  >] 


where 


y = I/I. 


TMs  equation  is  useful  for  evaluating  the  variance  of  log- 
irradiance  fluctuations  with  arbitrary  probability  distribution, 
because  the  histograms  of  irradiance  fluctuations  observed  often 
deviate  from  the  Idealized  log-normal  distribution.  The  variance 


of  log  irradiance  fluctuatons  for  the  case  in  which  the  Fresnel 
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zone  / XL  is  given  by 


c^(P)  = 3.2  1? 

2 

By  inserting  the  first-four  irradiance  moments  and  the  C = 
-11  -2/3 

2.98  X 10  m obtained  by  the  temperature  fluctuations 
measurements,  we  obtained  - 2.9  mm. 

Using  the  multi-path  optical  transmission  in  the  turbulence 
chamber,  we  have  measured  the  variance  of  log-irradiance  fluctua- 
tions  and  histograms  versus  optical  path  length  L,  The  prelimin- 
ary results  Indicated  clearly  the  saturation  effect  over  the 
longer  path  length.  This  was  presented  at  the  URSI  meeting  in 
Massachussets  (October  1976).^^  We  have  also  compared  the  direct 
transmission  measurements  of  irradiance  fluctuations  with  those 
using  the  flat  mirror  or  corner-cube  reflector.  According  to  the 
preliminary  results,  irradiance  fluctuations  of  collimated  He-Ne 
6328A  laser  beam  are  smallest  in  case  of  the  beam  reflected  by  a 
comer-cube  reflector  compared  to  the  direct  transmission  and  that 
reflected  by  a flat  mirror  over  the  same  optical  path  length.  The 
coefficients  of  variation  (Yq)»  skewness  and  excess  (Y2)» 

are  summarized  as  a table: 

corner-cube  reflector  direct  path  flat  mirror 


This  result  may  be  interpreted  as  follows:  Since  during  the  short 

transmission  time  the  turbulence  can  be  regarded  as  "frozen”,  the 
beam  reflected  by  a corner~cube  will  tend  to  go  back  following 
the  same  path.  Consequently,  forward  and  backward  transmission 
would  be  somewitat  correlated.  On  the  other  hand,  the  beam  reflected 
by  the  flat  mirror  would  propagate  in  much  wider  range  of  directions 
due  to  the  Snell's  law,  and  thus,  the  forward  and  backward  trans- 
missions are  statistically  more  independent  than  the  corner-cube 

reflector  case.  The  preliminary  results  were  presented  at  the 

* 12 

Optical  Society  of  America  Meeting  (October  1976).  We  plan  to 
analyze  these  experimental  results  more  precisely  by  formulating 
the  theoretical  model. 

Intensity  correlation  interferometry  using  the-  RETICON,  64 
linear  array  of  silicon  photodiodes,  (RL-64-P)  are  being  developed 
(cf.  Fig.  12)  in  order  to  clarify  the  Influence  of  atmospheric 
turbulence  on  the  intensity  correlation  interferometer.  Since  the 
size  of  the  diode  is  comparable  to  the  spacing  between  photodiodes, 
the  digital  signal  processing  for  reducing  the  directly  measured 
quantities  to  the  discrete  sampled  intensity  correlation  is 
essential.  A simple  case  was  analyzed  by  Hohamed  Ref ale,  a 
graduate  student. 

In  these  experiments,  the  influence  of  phase  fluctuations  due 
to  the  atmospheric  turbulence  is  an  important  subject  in  intensity 
correlation  interferometry.  Thus,  we  started  to  prepare  the  in- 
vestigation of  phase  fluctuations. 
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FIELD  EXPERIMENT 


A corner-cube  reflector  and  a flat  mirror  and  the  basis  for 

O 

a He-Ne  Laser  (Spectra  Physics  135^  6328A)  were  installed  on  the 
roof  of  the  Physical  Sciences  Building  (see  Fig.  13).  The  laser 
beam  from  Room  952,  Engineering  Building  was  transmitted  to  the 
reflectors  on  the  Physical  Sciences  Building  over  the  distance 
274  meter.  A small  receiving  telescope  (2"  diameter.  Zoom,  Swift) 
with  a photomultiplier  tube  (IP21)  with  an  interference  filter 
and  a pin  hole  mounted  on  the  micro-positioner  is  placed  next  to 
the  He-Ne  gas  laser  in  Room  952  (see  Fig.  14).  The  output  signal 
from  the  photomultiplier  tube  is  amplified  by  the  Tektronix 
differential  amplifier  and  monitored  by  a Cathode  Ray  Oscilloscope 
on  the  ninth  floor.  The  output  signal  from  the  differential 
amplifier  is  transmitted  to  the  Room  436  for  digital  signal 
processing  by  the  PDP  15/20  through  a 50  £i  coaxial  cable.  The 
wind  velocity  and  wind  direction  are  being  measured  on  the  roof 
of  the  building  and  their  information  is  being  monitored  at  the 
digital  data  acquisition  console.  The  details  of  digital  data 
acquisition  system  and  programs  used  are  identical  to  those  used 
for  laboratory  simulation  experiments  described  in  the  previous 
section. 

O 

When  a collimated  beam  (diameter  2”)  of  the  He-Ne  6328A  was 
transmitted  from  the  Engineering  Building  to  the  Physical  Science 
Building,  the  following  experiments  were  performed  in  the  evening 
of  14  Dec.,  1976  to  the  morning  of  15  Dec,,  1976  from  7:00  pm  until 
7:00  am.  The  reason  we  chose  night  time  was  that  tJje  atmospheric 
turbulence  is  relatively  steady  compared  to  the  day  time. 
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The  results  obtai.ied  are  compared  with  respect  to  coefficients 


of  variation  (Yq) , skewness  and  excess  (y^) J 


corner-cube  reflector 
0.71 


flat  mirror 
0.99 


direc^  path 
0.3826 


0.90 


0.66 


0.1239 


1.80 


0.79 


0.2045 


These  numbers  are  typical  examples  of  results  obtained.  These 
results  were  not  taken  simultaneously,  because  dedicated,  real- 
time data  processing  makes  simultaneous  experiments  impossible. 
During  the  digital  processing,  the  weather  conditions  might  have 
changed.  The  optical  path  length  for  the  direct  path  measurement 
is  the  half  of  the  total  path  length  of  either  comer-cube  or 
flat  mirror  reflector  cases. 

It  is  most  desirable  to  use  the  analog  instrximentation  tape 
recorder  for  storing  the  analog  data  of  direct,  corner-cube  and 
flat  mirror  cases,  simultaneotjsly  and  analyze  these  data,  later 
by  the  PDF  15/20,  For  this  purpose,  we  plan  to  utilize  the 
recorder  CEC  VR  3600,  (3  FM  Channels,  500  kc,  and  3 AM  Channels, 
1.5  MHz). 

We  have  installed  another  corner-cube  reflector  on  the  top 
of  the  campus  hill.  We  can  transmit  the  laser  beam  from  Room  952 
Engineering  to  the  reflector  over  the  distance  587.65  meters. 

We  have  taken  histograms  cf  irradiance  fluctuations  observed  under 
various  weather  conditions  and  also  measured  coefficients  of 
variation,  skewness  and  excess.  Comparing  these  results  with 
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log-normal,  Rlce-Nakagarai  distributions,  we  found  many  cases 

deviating  fro’i  these  distributions.  The  results  were  discussed 

12 

at  the  Optical  Society  of  America  Meeting  (1974).  This 
location  is  not  convenient  for  comparing  the  direct  laser 
transmission  because  of  the  lack  of  electric  power  supply. 
Therefore,  we  may  use  this  reflector  af«,^r  establishing  the 
model  of  corner-cube  reflector  in  conjunction  with  optical 
transmission  through  the  turbulent  atmosphere. 

During  the  year  of  1974,  we  installed  the  comer-cube 
reflectors  (6  units)  on  the  roof  of  the  former  Naval  Radar  site 
(N  - Calif  - 1015)  at  Pleasant  Peak  in  the  Cleveland  National 
Forest;  altitude,  4,000  ft.,  approximately  20  miles  on  the 
line  of  site  from  Room  436,  Engineering  Building.  Unfortunately, 
before  fully  utilizing  these  reflectors,  they  were  stollen.  Due 
to  inconvenience  and  lack  of  security,  we  have  discontinued  the 
use  of  this  site. 

During  the  year  of  1973,  we  obtained  the  coelstat  consisting 
of  the  16”  and  12  1/2"  flat  mirrors,  (Fig,  15).  The  16"  mirror 
can  be  controlled  remotely  so  that  we  may  seek  a specific  target. 
Wp  obtained  a 10"  Newtonian  telescope  matching  the  coelstat  and 
Installed  it  in  Room  952,  Engineering  Building  (Fig.  16).  These 
mirrors  and  telescope  were  manufactured  by  the  Cave  Optical  Co., 
Long  Beach,  CA.  After  completion  of  the  coelstat  telescope 
system,  we  transmitted  the  laser  beam  through  it  to  the  campus 
hill  and  analyzed  the  reflected  beam.  Then,  we  detected  the 
effect  of  the  minute  mechanical  vibration  due  to  the  disturbance 


of  wind  on  the  coelstat  by  analyzing  the  power  spectrum  of 
irradiance  fluctuations.  A good  wind  shield  such  as  prefabricated 


dome  is  necessary  to  utilize  the  system  most  meaningfully. 
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A CORNER  CUBE  AND  FLAT  MIRROR  REFLECTORS  AND 
A HE-NE  LASER  FOR  DIRECT  TRANSMISSION  EXPERIMENT 
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COELOSTAT  FOR  FIELD  EXPERIMENT 
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4.  COMPUTER  SIMUUTION 


In  order  to  facilitate  the  computer  simulation  of  optical 
propagation  in  the  turbulent  atmosphere,  we  have  developed  the 
following  two  basic  programs:  (1)  complex  Fourier  transform  and 
its  inverse  transformation  and  (2)  computer  generation  of 
gaussian  random  variables. 

The  diffraction  and  scattering  of  electromagnetic  waves  in 
turbulent  atmosphere  can  be  treated  most  conveniently  by  using 
the  complex  wave  amplitudes,  of  which  real  and  imaginary  parts 
are  mutually  Hilbert  transform.  This  method  is  well-known  in 
optics  as  the  analytic  signal  representation.  Using  the  PDP 
15/20  with  24  k words  coia  memory,  we  can  now  perform  compex 
Fourier  transform  of  512  words  of  complex  signals.  This  up- 
dated program  is  included  within  the  lECS  described  in  Appendix 
A,  This  program  was  developed  by  Wayne  Teeter,  a honor  student 
in  the  Senior  class,  as  an  engineering  aide  to  this  project. 

The  next  step  is  to  define  the  optical  transmission  function 
of  turbulent  atmosphere  and  to  calculate  the  diffracted  wave 
produced  by  a section  of  the  turbulent  atmosphere  at  the  observa- 
tion points.  The  numerical  analysis  of  diffraction  integrals 
can  be  handled  conveniently  by  successive  use  of  the  complex 
Fourier  transforms. 


The  optical  transmission  function  of  the  turbulent  atmosphere 
can  be  simulated  by  specifying  the  refractive  indexes  of  the 
turbulent  air  from  the  ensemble  of  computer  generated  random 
variables*  For  simplicity,  we  may  assume  the  gaussian  distribution 
for  refractive  index  at  least  at  the  beginning.  The  fluctuations 
of  the  index  of  refraction  should  be  non-gaussian,  because  the 
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refractive  Index  Is  proportional  to  the  temperature  of  which 
statistics  are  clearly  non-gauss Ian.  Thus,  we  have  been 
preparing  the  computer  generation  of  random  variables  which  can 
be  specified  by  the  central  moments  or  cumulants.  Since  we 
measured  up  to  4th  cumulants  in  the  laboratory,  computer  genera- 
tion of  random  variables  is  emphasized  on  the  first  four  cumu- 
lants . 

It  would  be  Interesting  to  note  that  the  complex  Fourier 
transform  will  also  be  useful  for  this  random  variable  genera- 
tlon  process.  For  Instance,  we  first  define  the  characteristic 
function  of  the  random  process  using  the  specified  cumulants 
and  then  calculate  the  probability  distribution  function  by 
applying  the  complex  Fourier  transform  to  the  characteristic 
function.  By  comparing  a computer  generated  random  nuinber 
X (0  - X - 1)  with  the  distribution  curve,  we  define  a 
random  variable.  This  is  essentially  an  encoding  scheme  to 
generate  signals  with  a non-gaussian  statistics  from  an  input 
pseudo-random  number. 

The  generation  of  gaussian  random  variables  using  the 
pseudo-random  numbers  generated  by  the  PDF  15/20  was  carried 
out  by  Loc  Ta,  graduate  research  assistant. 

It  is  an  interesting  and  useful  exercise  to  analyze  the 
case  of  corner-cube  reflector  and  the  plane  mirror  in  the 
propagation  through  the  turbulent  atmosphere.  We  may  establish 
models  for  corner-cube  and  flat  mirror  cases  by  comparing  the 
simulation  results  with  experimental  results  in  either  a 
laboratory  generated  turbulence  or  in  the  field. 
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CONCLUSION 


During  the  past  five  years,  we  have  studied  and  established 
various  basic  tools  for  systematic  investigation  of  Imaging  of 
incoherent  objects  through  the  turbulent  atmosphere;  pseudo- 
thermal  sources  for  incoherent  objects,  the  turbulence  chamber 
for  laboratory  simulation  and  computer  programs  for  statistical 
signal  processing,  automatic  measurements  and  generation  of 
random  variables,  A general  purpose  interactive  experimental 
control  system  program  has  been  developed  for  PDP  15/20  and  has 
been  fully  utilized  for  this  research  up  to  the  limiting 
capacity  of  the  computer.  These  accomplishments  will  be  useful 
and  essential  for  the  next  step  of  the  related  important  research 
such  as  establishing  the  models  of  the  turbulent  atmosphere  and 
feasibility  study  of  imaging  systems  for  incoherent  objects 
through  the  turbulence. 
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A GENERAL  PURPOSE  INTERACTIVE 
EXPERIMENT  CONTROL  SYSTEM  (lECS) 

The  lECS  is  based  on  the  Digital  Equipment  Corporation  PDP-15 
computer  with  2AK  words  of  main  memory,  two  magnetic  tape  drives 
(DEC-Tape),  a paper  tape  reader/punch  and  a Model  35  Teletype, 

In  addition  to  the  standard  hardware  several  "custom*'  items  have 
been  added.  They  are:  an  analog-to-digital  converter  (ADC)  with 

a sample-and-hold  front  end  and  a 7 microsecond  conversion  time, 

% 

two  digital  voltmeters  (DVM*s)  with  a 5 digit  BCD  readout  and  two 
digital- to-analog  (DAC)  converters  with  10-bit  resolution.  Fig. 
1,1  Illustrates  this  configuration  and  the  interaction  between 
the  devices.  In  support  of  the  hardware  attached  to  the  I/O  bus 
there  is  a graphics  pen-plotter  and  a CRT  display  oscilloscope 
attached  to  the  output  of  the  DAC's  which  provide  the  graphics 
output  for  the  lECS, 

The  control  of  these  add-on  peripherals  is  done  in  the  normal 
manner  through  input/output  commands  sent  along  the  I/O  bus. 
Several  machine  language  utility  programs  handle  the  transfer  of 
data  to  and  from  the  extra  peripherals  and  are  described  in  the 
following  section.  All  the  extra  peripherals  operate  under  the 
Interrupt  control  except  the  DAC’s  which  do  not  require  interrupt 
service.  The  relative  conversion  speeds,  of  the  input  peripherals 
(ADC  and  DVM’s)  is  such  that  an  overlap  in  the  sampling  speeds 
allows  a wide  range  of  signals  to  be  converted.  The  maximum 
sampling  rate  of  the  ADC  is  120KHz  while  the  conversion  rate  for 
DVM’s  is  approximately  10  to  200  milliseconds  depending  upon  the 
magnitude  of  the  volta^ic  being  sampled.  With  the  ADC  the  user 
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may  perform  most  sampling  tasks  currently  in  operation  in  our 
research  group.  The  sampling  speed  of  the  ADC  is  continuously 
variable  over  the  range  of  0.5Hz  to  20KHz  (or  fixed  at  27KHz  or 
120KH^  to  provide  for  the  variety  of  source  signals  present  in 
the  enviorment. 

The  graphics  output  devices,  pen~plotter  and  display  oscillo- 
scope, are  used  to  display  data  generated  by  direct'  acquisition 
or  indirect  computation  within  the  lECS.  The  graphics  pen-plotter 
and  the  display  ocsilloscope* are  controlled  using  the  10-bit  DAC's 
providing  a 1024  point  resolution  for  the  graphics  output  data. 
Whan  the  display  ocsilloscope  is  in  operation  the  current  image 
on  the  screen  is  refreshed  every  30’ th  of  a second  to  provide  a 
flicker-free  trace  for  viewing.  The  pen-plotter  in  turn  operates 
at  a much  slower  rate  and  is  capable  of  plotting  one  point  every 
80  milliseconds. 

The  standard  peripherals,  magnetic  tape  and  paper  tape  are 
serviced  through  the  normal  system  utilities  and  provide  a means 
for  mass  backup  storage  of  programs  and  data.  The  teletype  is 
used  chiefly  for  the  processing  of  language  commands  and  the 
printing  of  some  data  output,  A link  exists  to  the  PDP-10/70 
computer  through  the  magnetic  tape  where  any  large  printed  output 
may  be  taken  to  be  listed  on  the  high  speed  line  printer. 
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Given  the  best  hardware  system  available  a computer  still 
requires  software  to  implement  the  desired  actions.  It  is 
through  the  software  that  a computer  system  becomes  the  powerful 
computational  and  control  tool  that  we  witness  today.  The  soft- 
ware on  which  the  lECS  is  based  was  initially  developed  by  the 
Digital  Equipment  Corporation  for  use  on  the  PDP-8  computer.  The 
programming  language  FOCAL  was  designed  to  be  run  on-  the  first 
minicomputer  priced  in  the  range  affordable  by  the  general  public. 
This  computer  (the  PDP-81)  had  a memory  of  4K  words » a Teletype 
and  a paper  tape  reader/punch.  It  was  a "bare  bones"  system  to 
say  the  least.  To  be  of  any  use  to  our  application,  the  FOCAL 
language  had  to  be  extended  to  include  control  of  the  peripheral 
devices  attached  to  the  I/O  bus.  The  original  PDP-15  version  of 
the  FOCAL  interpreter  was  first  extended  by  several  individual 
groups  including  the  University  of  Mississippi  Medical  School 
and  a group  of  users  at  the  Los  Alamos  Scientific  Laboratory. 

The  current  version  of  FOCAL  operating  in  conjunction  with  the 
special  utilities  package  that  comprises  the  lECS  is  a modified 
version  of  the  above  two  enhanced  versions.  It  would  be  correct 
to  say  that  the  FOCAL  processor  operating  under  the  lECS  is  a 
far  cry  from  the  original  FOCAL  language  processor  first  imple- 
mented on  the  PDP-8  in  1965. 

Fig,  1,2  outlines  the  basic  software  configuration  as  cu  "ent- 
ly  Implemented  at  the  School  of  Engineering  of  the  University  of 
California,  Irvine,  The  system  consists  of  a FOCAL  language 
interpreter  which  is  used  to  process  the  FOCAL  language  programs 
generated  by  the  user,  edit  these  programs  using  the  Teletype  and 
provide  an  overall  control  over  the  operational  structure  of  the 


lECS.  Working  In  conjunction  with  the  FOCAL  interpreter  is  the 
utilities  pavkage  which  transforms  the  original  FOCAL  interpreter 
into  the  lECS.  This  package  consists  of’  special  functions  needed 
to  provide  full  control  over  the  data  acquisition  and  display 
processes  (the  computational  ability  is  contained  within  the 
original  interpreter).  The  third  major  portion  of  the  lECS  is 
the  run-time  software  required  to  interface  to  the  mass  storage 
devices,  control  the  Teletype  and  provide  floating  point  com- 
putation capabilities  needed  during  the  execution  of  FOCAL 
language  programs  (the  current  installation  of  the  PDP-15  does 
not  have  a hardware  floating  point  processor) . 


HARDWARE  CONFIGURATION 


PDP  15/20  DATA  ACQUISITION  SYSTEM 


REVIEW  OF  PERTINENT  RESEARCH 


A pseudo-thermal  light  source  was  used  for  simulating  the  Brown- 1V?iss 

intensity  interferometer  by  Martienssen  and  Spiller^  and  later  by  Haner  and 
2 

Isenor.  From  the  intensity  interferogram  in  the  form  of  the  transparency 

obtained  by  the  system  similar  to  those  mentioned  above,  the  image  of  a 

pseudo-thermal  light  source  was  reconstructed  by  applying  the  Fraunhofer 

3 

diffraction,  an  optical  analog  of  Fourier  transform.  The  design  and 
results  of  Narrabri  intensity  interferometer  constructed  in  Australia  by 

4 

R.  Hanbury  Brown  was  described  in  his  book  on  the  intensity  interferometer. 

In  this  book,  Hanbury  Brown  emphasized  the  freedom  of  the  system  from  the 
atmospheric  scintillations  as  a remarkable  and  important  property  of  the 
interferometer  and  mentioned  that,  *’It  would  be  worthwhile  to  refine  both 
the  analysis  and  the  observations  to  establish  these  effects  more  precisely." 
(P,  72,  Reference  A). 

The  unique  determination  of  the  image  of  an  arbitrary  incoherent  object 
is  not  possible  by  using  the  Hanbury  Brown  intensity  interferometer  because 
the  phase  of  mutual  intensity  cannot  be  obtained.  One  of  the  interferometric 
systems  which  can  recover  the  phase  information  is  Jennison's  phase  sensitive 
interferometer,  which  measures  interferences  between  three  beams  of  incident 
radiation  field. ^ Both  R.C.  Jennison  interferometer  and  the  triple  intensity 
are  not  sensitive  to  the  atmospheric  phase  disturbances.  Further  detailed 
Investigation  of  this  effect  is  worthwh-* Te. 

Another  type  of  intensity  interferometer  which  will  not  be  disturbed  by 
the  atmospheric  phase  fluctuations  and  can  uniquely  reconstruct  the  incoherent 
object  is  the  intensity  interferometer  with  a coherent  background  due  to  the 
principal  investigator.  The  astronomical  interferometer  using  the  heterodyne 
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detection  is  intimately  related  to  this  interferometer.  Tlie  heterodyne 
Interferometer  for  the  astronomical  observation  was  studied  by  Nieuwen- 
huijzen,^  Van  de  Stadt^*®  and  Johnson,  et  al.^  The  laboratory  experiment 
also  related  to  the  interferometer  was  performed  by  Soemarjona  and  Hirano,^^ 
The  compound  intensity  interferometer  which  measures  the  intensity  corre- 
lation between  one  beam  and  the  product  of  two  beams  with  a small  separa- 
tion proposed  by  Me  Phie^^  is  also  pertinent  to  the  subject  of  this 
research  grant. 

The  magnitude  of  mutual  intensity  due  to  an  incoherent  source  such 

as  a star  can  also  be  determined  by  measuring  the  ensemble  average  of 

spatial  autocorrelation  functions  of  the  instantaneous  image  intensity 

under  the  influence  of  atmospheric  turbulence.  This  was  pointed  out 

12  13 

independently  by  Labeyrie  and  Asakura,  et  al.  This  method  is  called 
the  speckle  interferometry  and  has  recently  been  reviewed  by  Dainty. 

Note  that  the  phase  Information  of  mutual  intensity  cannot  be  determined 
by  the  speckle  interferometry  in  general, 

• Certain  classes  of  pattern  of  incoherent  light  sources  can  be  deter- 
mined by  observing  only  magnitude  of  mutual  intensity  by  the  Hanbury  Brown 
and  Twiss  intensity  interferometer,  if  the  phase  of  mutual  intensity  can 
be  derived  from  its  magnitude  uniquely  or  phase  information  of  certain 
patterns  are  predetermined.  The  relations  between  phase  and  magnitude 
of  the  mutual  intensity  are  given  by  the  Hilbert  transforms.  Phase  problems 
in  coherence  and  theory  and  Intensity  interferometry  were  studied  by 
Nussenzvieg^^  and  Bates. It  should  be  noted  that  some  results  on  the 
analogous  relations  extensively  studied  in  the  electric  network  theory  may 
become  very  useful  in  optical  problems.  Some  symmetric  patterns  which  can 
be  detected  without  additional  phase  information  were  studied  by  the 
principal  investigator  in  conjunction  with  the  optical  pattern  recognition.^^ 
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The  image  reconstruction  from  the  modules  of  mutual  intensity  function  were 

18 

studied  by  Kohler  and  Mandel  using  the  analytical  properties  of  mutual 

19 

intensity  function.  The  phase  determination  of  mutual  intensity  in  the 
presence  of  atmospheric  seeing  has  been  discussed  by  Rogstad  with  respect 
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to  multi-element  interferometry  along  the  line  of  Jennison  interferometry. 

21 

An  optical  telescope  array  for  intensity  interferometer  was  treated.  A 

method  for  recovering  phase  information  which  is  missing  in  the  speckle 

interferometry  was  treated  by  using  the  autocorrelation  of  Fourier 

transform  of  image  intensity  influenced  by  the  atmospheric  fluctuations. 

The  phase  errors  in  Michelson  interferometer,  their  indentification  and 

the  removal  was  treated  by  using  the  concept  of  complex  zeros  of  the 
23 

interferogram.  A more  unified  approach  to  the  imaging  viewed  through 

randomly  fluctuating  media  based  on  interferometry  and  holography  was 

..  j 24 
presented. 

Important  review  articles  and  books  on  wave  propagation  in  the 

turbulent  atmosphere  are  referenced  in  24-29.  The  English  translation 

25 

of  V.I.  Tatarski  report  has  been  utilized  extensively  in  our  recent 
research.  Reference  27  is  useful  to  provide  the  up-to-date  overvie\'?  of 

current  research  activities  on  the  optical  propagation  in  the  atmosphere 

29  2 

in  many  countries.  The  forthcoming  book  by  A.  Ishimaru  includes 
in  the  turbulent  atmosphere  due  to  the  solar  wind  as  well  as  the  most 
up-to-date  review  of  multiple  scattering.  The  wave  propagation  in  the 

turbulent  atmosphere  plays  an  essential  part  in  adaptive  optical  system, 

30 

which  is  very  well  covered  in  the  recent  special  issue. 
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